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If massive neutrinos possess magnetic moments, a magnetic field 
can cause a spin flip. In the case of Dirac neutrinos the spin flip 
converts an active neutrino into a sterile one and vice versa. By 
contrast, if neutrinos are Majorana particles, a spin flip converts 
them to a neutrino of a different flavor. We examine the behavior 
of neutrinos in a random magnetic field as it occurs, for instance, 
in certain astronomical objects, such as an active galactic nu- 
cleus. Both Dirac and Majorana neutrinos behave ergodically: 
independently of their initial density matrix, they tend towards 
an equipartition of the helicity states. As a result, about half 
of the Dirac neutrinos produced becomes sterile. For Majorana 
neutrinos, there will be an approximate equipartition of flavors, 
independently of the production mechanism. 
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According to the standard model of electroweak interactions, neutrinos can- 
not possess magnetic moments. The possibility of giving magnetic moments 
to neutrinos arises in minimal extensions of the standard model: once right 
handed neutrinos are introduced, chirality violating, Pauli type magnetic 
interactions with an electromagnetic field can be introduced. There is an 
important difference between the physical effects caused by such interac- 
tions, depending on whether neutrinos are Dirac or Majorana particles. The 
Dirac case is somewhat uninteresting: a chirality flipping interaction causes 
a conversion between active and sterile neutrino species. By contrast, in the 
case of Majorana neutrinos, different chiralities must correspond to different 
flavors and thus, a chirality flipping transition causes transitions between 
flavors at the same time. The minimally extended standard model gives rise 
- via higher order loop effects - to anomalous magnetic moments; however, 
their magnitudes are rather small and one can hardly hope to ever detect 
them. One needs truly "new physics" beyond the standard model in order 
to generate detectable magnetic moments of neutrinos. There exist rather 
strong - albeit somewhat model dependent - upper limits of flavor diagonal 
magnetic moments, typically in the range of 10~ 12 /iB, where \xb is the Bohr 
magneton. The upper limits for transition moments are generally weaker, 
typically by some three to four orders of magnitude; see, for instance, ref. 
for a recent review. Thus, the experimental identification of any effect caused 
by the presence of a magnetic interaction of neutrinos is of utmost interest, 
for it provides a handle on physics beyond the standard model. 

In this work we examine the fate of neutrinos possessing anomalous mag- 
netic moments in random magnetic fields. Such fields may serve as models 
of the chaotic fields expected to exist within certain astrophysical objects, 
in particular, within active galactic nuclei (AGN). An AGN is a particularly 
interesting environment, since, according to many speculations, it is a source 
of ultra high energy neutrinos, see, for instance the proceedings of the Hawaii 
meeting 0. 

It is generally assumed that most neutrinos emerge from the production 
and subsequent decay of hadrons. Assuming that one can extrapolate accel- 
erator based data to the relevant (PeV or higher) energies, the overwhelming 
majority of hadrons produced in a hadron collision consists of pions. Hence, 
the charged leptons produced in a neutrino telescope observing the point 
source must be mostly muonsQ. We argue that, under appropriate circum- 
stances, a random magnetic field equilibrates the helicities; hence, in he case 

2 In a previous work || we discussed the possibility that resonant interactions in the 
hot plasma present in an AGN may cause the appearance of different flavors of neutrinos. 
Also, Learned and Pakvasa M considered flavor conversion due to long range oscillations. 
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of Majorana neutrinos, the flavors as well. In this work, we consider the 
behavior of a single spin-1/2 field, without paying detailed attention to the 
flavor structure. 

In order to describe the average behavior of a neutrino in a random mag- 
netic field, one has to solve the dynamical equations governing the propa- 
gation in an arbitrary magnetic field. The solution then has to be averaged 
over the ensemble of magnetic fields. 

We use the front form of dynamics |J . As explained in a previous work || , 
this formulation of dynamics is advantageous in a situation in which one 
considers the propagation of high energy particles (E 3> m, where m is the 
rest mass) and in which certain discrete symmetries, such as C and P play 
no significant role. Clearly, the propagation of high energy neutrinos falls 
into this category. 

We begin with the usual Dirac Lagrangian of a particle in an external 
electromagnetic field, F^ v : 

L = i> (i<fd» + m + ^F""<v) V (1) 

We work in the rest frame of the magnetic field. Assuming the field to be 
a static oneQ, we can set F 0i = 0, = e^B^. In the case of interest one 
has to solve the Dirac equation in an arbitrary static magnetic field, since we 
want to average the solution over an ensemble of the No explicit solution 
is known for such a problem. However, we proceed to show that in the high 
energy limit the problem can be solved in a closed form. 

We introduce a coordinate system in which two of the coordinates are null 
directions corresponding to characteristic lines of a relativistic wave equation, 
viz. : 

t = — (x° - x 3 ) , z = -= (x° + x 3 ) and x A ; (A =1,2). (2) 

Correspondingly, the metric is of the form, 

9zt = 9tz = l, 9ab = Sab, (3) 

and all other components vanish. 

A Dirac spinor can be decomposed along the null directions given in (||]) 
by introducing the mutually orthogonal projectors, 

Pt = \ltl\ P z = lzl z (4) 



3 In physical terms, this means that the characteristic time scale of change of the field 
is large compared to the time of passage of neutrinos. 
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In what follows, we use the shorthand, 



(5) 



It is a straightforward matter to decompose ([!]) according to the conju- 
gate null directions and express it in terms of the variables <fi and x- The 
purpose of such an exercise is a very simple one. If, for the sake of definite- 
ness, t is regarded the "time" variable describing the dynamics of the system, 
only <p obeys an equation containing d t . Hence, the component of the Dirac 
spinor corresponding to the conjugate null direction obeys only an equation 
of constraint. The constraint can be, in turn, solved before one attempts to 
attack the problem of dynamics. 

After carrying out the decomposition of (|IJ) according to the null direc- 
tions, one finds: 



L = V2 [<jf (id t - iV2^e AB lA B B 
+ X ] (id z -iV2iie AB lA B B )x 



+ 



V2 



V ( h A d A +m- ^=l*B 3 t A Bl A l B \ X 



+ X V Ul A dA + m- e ABl A -/ B 



(6) 



Variation of (|6|) with respect to x^ gives the constraint. The constraint 
can be solved in a straight forward fashion and eliminated from the La- 
grangian. The result is conveniently written in Hamiltonian form: 



L 
H 



+ 



-ia B e BC p c + m — jj\/2B 3 a 3 



x Q -B 



x 



(io R e RS p s + m- nV2B 3 a3 



(7) 



Solving the constraint eliminates two components of the original, four compo- 
nent Dirac spinor. Therefore, instead of the original Dirac matrices one can 
use 2x2 Pauli matrices. One easily verifies that —ie AB j B — > a A gives the cor- 
rect representation. We also introduced the Hermitean operator, p A = —id A 
for the transverse degrees of freedom. 

The canonical momentum is given by n = i\pl<$ . (Of course, the odd 
looking factor of \/2 in the definition of the canonical momentum can be 
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eliminated by rescaling the time variable.) In equation (^), Q is an operator 
with matrix elements: 



(z\Q(B A )\z') =^=exp^V2£dz'e ABl A B B ^e(z-z') (8) 

All symbols of integration over z have been omitted. Eq. (^) is local in t 
and x A ; those arguments have been suppressed. 

The Hamiltonian appearing in (|7|) is exact. However, it is given by a 
rather complicated, non local and non linear expression: this is the cost we 
have to pay for explicitly eliminating the constraint. We now argue that 
one can introduce physically reasonable simplifications, as a result of which 
the problem becomes a manageable one. First of all, we notice that the 
exponential appearing in (|8]) is of modulus one. Furthermore, it is recognized 
as the matrix element of an operator which is unitary and Hermitean at the 
same time, hence its eigenvalues are ±1. Hence, one expects that at large 
values of \z — z'\ the exponent oscillates rapidly and thus contributes little 
to the Hamiltonian. The dominant contribution is thus coming from small 
values of the difference of longitudinal coordinates. Hence, it is reasonable 
to approximate the exponential in (|5[) by 1. In the remaining expression, 
one term is local in all variables and the remaining ones are proportional to 
e (z — z'). Hence, in a Fourier representation, viz. upon writing 

(J) (t, z, x A ^J = J dkip (t, k, x A ^ exp (— ikz) (9) 

and 

£(2) = l7/f exp( - lfc) - (10) 

(In the last equation V stands for the principal value.) Hence, at high energies 
(k ^> m) the Hamiltonian can be approximated by the local term. 

Neglecting terms of O (A; -1 ), the equation of motion for the density matrix 
in coordinate representation reads: 

— id t (z,x\p(t) z',x') = pV2a ■ B ^x, -j= \ ( z, x\p{t) z', x') 

- pV2(z,x\p(t)\z',x')aB [x',^J^ (11) 

In this equation, x stands for the transverse part of the coordinate and 
a ■ B is the two dimensional scalar product in transverse space. Of course, 
the coordinate x 3 had to be expressed by z and t; hence the t-dependence in 
the magnetic field. 
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We choose the initial condition so as to describe a neutrino produced at 
x = and with a fixed value of k: 



(z,x\p(0) z'x') = 5 2 { 



x) 5 2 (x 1 ) 



exp ik (z — z') 

2nk 



Ps (0) , 



(12) 



where p s (0) is the initial value of the spin density matrix. 

The variable k being large, the function exp ik (z — z') is rapidly oscillat- 
ing unless z ~ z' . Therefore, it is permissible to put z — z' in the coefficient 
of the exponential in (0). Further, in the approximation used, the dynamics 
described by eq. flTT| ) is independent of k and of x. Therefore, the dependence 
of p(t) on k and x is entirely determined by the initial condition. Thus, the 
dynamical equation reduces to an equation involving the spin density matrix 
alone, as in non relativistic spin dynamics. Thus, from now on, we omit the 
the subscript s and we have: 



— id t p (t) = \x\Pl 



a 



B 



z-r 



pit) 



(13) 



(Here and in what follows, x = is understood.) 

This equation can be solved by the standard time ordered series, viz. 



p(t) = p(0) 

+ ipV2 { dt' 
Jo 



a-B 



'z-r 

. V2 , 



,p(o) 



+ 
+ 



2! 



dt'dfT 



a-B 



'z-r 

. V2 , 



a-B 



'z-r 

. V2 , 



,P(0) 



(14) 



We choose the initial condition as: 



p(0) = -(l + Sa 3 ), 



(15) 



since neutrinos are produced with a definite helicity. (In the case of Dirac 
neutrinos, S — ±1, depending on whether a neutrino or anti neutrino is pro- 
duced. In the case of Majorana neutrinos, 5* may assume any value between 
the limits stated above, depending on the production mechanism.) 

Next, we average the solution, (|14"1) over the magnetic field. We choose 
the generating functional of the moments as follows: 

Z[j] = JvBexp- - J d^xd^x'B, (x) (x, x') B 3 [x') 
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x exp J d 3 xji (x) Bi (x) 



C 



-i 



L 



An(B 2 } 



Sij 



V 2 . 



L 



-2 



V 



2\ 2 S 3 



x — x') . 



(16) 



In the last equation, L and (B 2 ) stand for the correlation length and mean 
square magnetic field, respectively. The measure is normalized such that 
Z [0] = 1. The transverse projector is needed in order to make the correlation 
functions solenoidal. With the choice of the tensor C _1 given in (|16D , the 
leading term in the long distance behavior of the correlation function is oc 
exp — \x — x'\. In order to average equation (|T4]) over the magnetic field, one 
integrates over B 3 and sets the third component of the source equal to zero. 
The transverse generating functional reads: 



VB exp 



L 



8n(B 2 } 



d x 



B A (x) 8ar - 



d A d B ( x) 

^2 

V 



B 1 



x exp i / d 3 xj (x) ■ B (x) 



(17) 



We now notice that in the equation fll4|) , terms containing odd powers of 
\i are also odd in B A . Therefore, in the limit j —> the average of those terms 
vanishes. The even terms in the series are obtained by taking the appropriate 
functional derivatives of (p]). All of them are expressed in terms of multiple 
time integrals of (\t — 1'\) and its powers: those integrations are easily 
performed. It is sufficient to illustrate the procedure for the second order 
term in (|14]). 

Carrying out the integrations, one gets: 



-ti 2 -S( f dt'dt" 
2 Jo 



a 



B, 



a- B,a 3 ) = -ii 2 a 3 (B 2 )tL 1 - exp 



For large times the result in the last equation is just proportional to t. The 
higher order terms follow a similar pattern. The end result is: 

(p(t)) ~ g ( X + 5a 3exp-^ 



with 



^ = 2^ 2 (B 2 )L. 



Thus we arrive at the remarkable result that in the random field the 
behavior of the helicities is an ergodic one: irrespective of what the initial 
density matrix was, for t ^> T, the helicities are equally distributed. In the 
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case of Dirac neutrinos, this is rather uninteresting: roughly 1/2 of them 
is a sterile one. However, calculated neutrino fluxes emerging from such 
astrophysical objects as an AGN usually cannot be trusted to an accuracy 
which would permit an observational testing of the result. The presence of 
sterile Dirac neutrinos does, however, play a role in the early universe, see ; 
this paper contains a virtually complete bibliography on the subject^. 

By contrast, in the case of Majorana neutrinos, the two helicity states cor- 
respond to two different flavors. Given the fact that the neutrinos produced 
arise mostly from pion decay, the presence of the other flavor in roughly 
equal proportion is an observationally testable result. We conjecture that 
the situation is similar if all neutrino flavors are properly taken into account. 

This work was done during the authors' visit at the Dipartimento di 
Fisica, Universita di Firenze. We wish to thank Roberto Casalbuoni, Director 
of the Department for the hospitality extended to us. We also thank Bianca 
Monteleoni for useful conversations on observational neutrino astrophysics 
and to Kari Enqvist for some critical remarks. 
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